Dynamical Density Fluctuation of Superfluids near Critical Velocities 



O 

1—5 

in 



I 

o 



> 

0^ 
0^ 

O 
O 



X 



Yusuke Kato^ and Shohei Watabe^'B 

^Department of Basic Science, The University of Tokyo, Tokyo, 153-8902, Japan 
^Institute of Physics, Department of Physics, The University of Tokyo, Tokyo 153-8902, Japan 

We propose a stability criterion of superfluids in condensed Bose-Einstein systems, which incor- 
porates the spectral function or the autocorrelation function of the local density. Within the Gross- 
Pitaevskii-Bogoliubov theory, we demonstrate the validity of our criterion for the soliton-emission 
instability, with use of explicit forms of zero modes of the Bogoliubov equation and a dynamical 
scaling near the saddle-node bifurcation. We also show that the criterion is applicable to the Landau 
phonon instability and the Landau roton instability within the single-mode approximation. 
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Introduction. Studies of the critical velocity, above 
which a superfluid becomes unstable, can trace its his- 
tory back to the discovery of superfluidity in ^He 
The critical velocity is determined by either of the fol- 
lowing, whichever yields the smaller critical velocity: the 
instability of the excitation spectrum with respect to a 
flowing superfluid (the so-called Landau instability Q) 
or the emission of topological defects such as quantized 
vortices [l,!!]. 

The realization of Bose-Einstein condensation in cold 
atoms has renewed interest in this issue. The critical 
velocities in cold atoms have been measured by moving 
a blue-detuned laser beam in a Bose-Einstein conden- 
sate (BEC) where superfluidity was found to break 
down with the emission of vortices at the velocities much 
smaller than the Landau critical values. Those experi- 
mental results Q are consistent with numerical results 
of the time-dependent Gross-Pitaevskii (GP) equa- 
tions [13, El- 

From the viewpoint of non-linear physics, the instabil- 
ity of superfluidity found in lid , [llj was identifled as a 
saddle-node bifurcation 12, 13|, where a stable and an 



unstable steady solutions for the GP equation merge [14 1 . 
As a simpler analog of the vortex-emission instability, 
the instability of the one-dimensional superflow against 
a penetrable potential barrier has been intensively stud- 
ied [l5-17| within the GP equation; superfluidity breaks 
down with the emission of solitons much below the Lan- 



dau critical velocity at a saddle- node bifurcation [15|, [18| . 

If we could find a common aspect between the two 
different kinds of instability: the Landau instability and 
the saddle-node bifurcation, it would provide a crucial 
step toward a coherent understanding of various kinds of 
instability in superfluids. In this Letter, we focus on dy- 
namical density fluctuations (DDEs) and show that the 
DDEs of superfluids are enhanced near the critical veloc- 
ity both in the Landau instability and soliton emission 
instability. On the basis of this finding, we propose a 
criterion for stability of superfluids fio']. 

Model. We quantify DDEs through the spectral func- 
tion of the local density, since the dynamical structure 
factor, usually denoted as S'(q, w) in the literature, is 



not appropriate to measure density fluctuations in the 
presence of obstacles. The spectral function of the local 
density is given by 

/(r,a.) = \{M\^)H^)\g)?5{^ -Ei+ E,). (1) 



Here, |g) is the state vector of the ground state or a stable 
superflow state, the energy of which is denoted by Eg. Ei 
denotes the energy of an excited state I. The summation 
runs over the excited states. 

We consider Bose systems with a short-range repul- 
sive two-body interaction at zero temperature. The x- 
axis is taken to be the direction of the supercurrent. 
We introduce a penetrable and short-ranged potential 
barrier Ucxix), which is localized along the x-direction 
around x = (independent of y and z). The shape 
of the barrier is similar to that shown in Eig. 1 a) of 
[20| . We treat the condensate and excitations within the 
scheme of the GP [1, @| and the Bogoliubov 21 1 equa- 
tions, respectively. The field operator is thus given by 



Tp{r,t) = *(r) + Ej-KWaj-e-"^^* - v* {r)a]e'^^*]. The 
function ^(r) represents the wavefunction of the conden- 
sate and satisfies the stationary GP equation: 



+ t/e.t(x)-M+l*(r)P 



*(r)=0. (2) 



The symbols aj and denote, respectively, the annihila- 
tion and creation operators of a one- particle state j. The 
wavefunctions Uj{r) and Vj{r) of the state j satisfy the 
Bogoliubov equation: 
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+ C/cxt(a;)-M + 2|*(r)|2 



(3) 
(4) 



and also satisfy the ortho-normal condi- 
tion / dr[uj(r)u*.(r) - Uj(r)u*(r)] = 5jk, and 
J dr[uj(r)ufc(r) — Wj(r)ufe(r)] = 0. We have used 
dimensionless units in ^ and 



2 



In the presence of a potential barrier Uck{x), we con- 
sider ([2]) under the boundary condition 



In terms of M{x, to), the spectral function is expressed as 



±00) 



exp 



i(Vx + ip±^ 



(5) 



The quantity ip = ip-^-—ip^ represents the phase difference 
between the two condensates separated by the potential 
barrier. V represents the superfluid velocity far away 
from the potential barrier. The relation between V and 
9? corresponds to the Josephson relati on [ig . |22| and its 
typical profile is given, e.g. in Fig. 2 of [16|. For chemical 
potential /i, the relation /i = 1 + follows from the 

boundary condition ([5]). 

When Ucx{x) — Uo6{x) with Uo > 0, the following 



facts are known from the analysis of Refs. [15|, [171, [1 



(a) . The critical velocity Vc is smaller than the Landau 

critical velocity(=l). 

(b) . When < V < Vc, a. stable and an unstable steady 

solutions to ([2]) exist. 

(c) . When V Vc — 0, the stable and the unstable 

steady solutions to ([2|) merge. 

As for V > Vc, no steady solutions to ^ with dissipa- 
tionless flows exist. When V is slightly above Vc, the 
emission of gray solitons is observed in numerical calcu- 
lations of the time-dependent Gross-Pitaevskii equation 
15|. When Ucxix) is a rectangular barrier, a Josephson- 
type V-(p relation was also derived in |22|. A rectan- 
gular barrier describes qualitatively generic properties of 
a penetrable short-range barrier. We thus assume the 
V-if relation for a generic short-range barrier to be a 
Josephson- type relation, from which (a), (b) and (c) fol- 
low. 

DDF in the presence of a potential barrier. The spec- 
tral function I{x, lu) at low cj is derived via the method 
used in (23|; first we find the zero modes of the Bogoli- 
ubov equation ([3]) and next obtain the solutions of ^ in 
the form of the series with respect to energy. On the ba- 
sis of the low energy asymptotics of the matrix element 
and the density of excited states, we discuss the details 
of I{x,uj) at low uj. 

In the Bogoliubov theory, the spectral function ([T|) re- 
duces to/(r,a;) = J2j 

For simplicity, we first discuss the case d = 1. When the 
volume fl of the system is sufficiently large, a typical 
value Ae of the level-spacing 
ri~^ for small e,-. 



Ej I is of the order of 



We here introduce a function M{x,uj) satisfying 
M{x,uj) = fl\-^*{x)uj{x) - ^{x)vj{x)\^, 



(6) 



when UJ — Sj for each j. Further we require that AI{x, uj) 
varies over the energy scale much larger than Ae, i.e. 



dM{x,uj)/duj Ae < M{x,uj) . 



(7) 



I{x,uj) = M{x,uj)D.'' 



i 



5{(jj - Sj) = M{x,uj)D{uj). 



(8) 

Here D{uj) denotes the one-particle density of states for 
the Bogoliubov mode, which becomes a constant at low 
UJ for c? = 1 in the thermodynamic limit {fl 00). 

In order to discuss M(x,uj) at low uj, we first recall 
that {u{x),v{x))*' = = -00 is a solution 

of ([3]) with zero eigenvalue (zero mode) for V < VC.[2^. 
Next, one can verify from the following observation that 

{u{x),v{x)Y = {d'f{x)/dip,-d'f*{x)/d(f)'=i()c (9) 



IS an additional zero mode of © at F = K[il|2i|. Re- 
garding y as a function of ip and taking the derivative of 
([2| with respect to ip, we obtain 

Kid'i'/dif) + ^^{d-i'*/d^) = V{dV/d^)-f. (10) 

This equation and its complex conjugation are put to- 
gether as 



C^c = V{dV/d^)'^o- 



(11) 



When V — Vc, V is maximal and hence dV/d(p = 0. The 
right-hand side in (fTTj) then vanishes and the proof is 
completed. 

When V < Vc, the wave function for a state j with a 



small e," can be shown to have the form 23l 2 



{uj{x),Vj{x)Y = (ejil) i {cipoix) + cEjipiix) + 0{e^j)} 

(12) 

Here ipi (x) is the solution of £Vi — ■00 which does not 
diverge exponentially at large \x \. When V 
other hand, we have [l 



23, 2 



Vc, on the 



{uj{x),Vj{x))' = {e^n) 5 {cipoix) + c'lpcix) + 0{ej)} . 

(13) 

Here c and c' are constants. 

tpo represents a phase fluctuation and hence does 
not contribute to density fluctuations, while Vc does . 
Correspondingly, the matrix element M{x, uj) behaves as 



ujf{x) 



M{x,uj) M{a 



0|'^(z,y)r 
dip 



V<Vc. 



V = Vc. 



(14) 

at low UJ. Here f{x) represents the contribution from ipi. 
The low UJ behavior of I{x,uj) for d = 1 can be deduced 
as M{x,uj)D{u}) from ^ and ([T4)) . 

The generalization to d = 2 and 3 is straightfor- 
ward as outlined below |26|; we classify the set of the 
eigenstate j with the angle (j) between the wave vec- 
tor and the potential wall, and introduce M{x,uj,4>) 
satisfying ^ and ([7]) for j with an angle 4>. At low 



3 



w, we can show that M{x,uj,4>) becomes independent 
of (f) and that it approaches M{x,u!). Consequently, 
I{x, Lu) — > M{x, uj)D{uj) at low w. Since D{uj) oc w'^'^^ at 
low uj for the Bogoliubov mode in d{= 1, 2, 3)-dimcnsion, 
we obtain 



d-2 



dip 



V = Vr. 



(15) 



at low UJ. Here fd{x,V) is a function of x and I'd is a 
constant. 

In (jl5p . we note that the exponent in cj jumps from d 
to d—2 when V approaches Vc from below. This behavior 
can be explained by a dynamical scaling relation near the 
saddle-node bifurcation fl^, Near the saddle- node 
bifurcation, there is a characteristic frequency a;*, which 
scales as w* oc \V-Vc\^. In [ill, it 

was shown that the 
decay rate F of an unstable steady state slightly below 
Vc and the emission rate F' of gray soliton slightly above 
Vc follow the scaling relations equivalent to F, F' oc \V ~ 
14 1 5 [g^]. We thus anticipate that uj* yields a crossover 
behavior in the w-dependence of I{x,uj) near Vc- 
Accordingly, we put I(x, cj) in the scaling form: 



I{x,u) =Lu''-'F{x,Lo\V -Vc\—^), 

where the asymptotic behaviors oi F{x,uj = uj\V - 
are determined from (1151) as 



F{x,uj) 



Vd 



fd{x) 

dip 



for cZi 0, 

for a) —> cx), 



(16) 



(17) 



with fd{x) = (Vc - V)fdix, V). The scaling relation ^ 
with (ITT)) shows that I{x, uj) is enhanced near V = Vc at 
low UJ. We have confirmed this scaling relation numer- 
ically. Figure [1] shows the scaling plot uj'^~'^I{x = 0,a;) 
vs. uj for Uex{x) = UoS{x) with Uo = 10 near Vc. We see 
that for each dimension, the data for various values of V 
near Vc collapse onto a single curve. 

In deriving the scaling relation ([16]) with (|17p , the fol- 
lowing facts were crucial: (I) the zero mode which cou- 
ples to density fluctuations exits only at V = Vc, and 
(II) the characteristic frequency uj* scales as \Vc — V\^^^ 
near the saddle-node bifurcation. These two facts hold in 
superfluid Bose systems where the critical velocity is de- 
termined by a saddle-node bifurcation. We thus expect 
that a scaling relation similar to (fTBl) holds also in the 



Hill], where the 



case of a superfiow around a disk [IC 
vortex-emission instability of superfluids has been iden- 
tified with a saddle-node bifurcation. 

DDF near the Landau critical velocity. We next show 
that the enhancement of the DDF near V — Vc also oc- 
curs in the Landau instability. To that end, we consider 
the DDF near the Landau critical velocity in spatially 
uniform systems within the single mode approximation 




FIG. 1: (Color online) Scaling relation near the saddle-node 
bifurcation for dimensions d = 1, 2 and 3. The scaled spectral 
functions F{x = 0,Li = - Kl"^/^) = ij^-''l{x = 0,uj) 
near the critical current are shown as a function of scaled 
frequency iD. The critical velocity Vc in dimensionless form is 
given by K = 0.049753 • • • for Ue^x) = UoSix) with Uo = 10. 
Circle (red), square (blue), and triangle (green) points are 
data for V = 0.04975, 0.0497 and 0.049, respectively. 



(SMA) |2S]. We assume that the excitation energy Sq in 
the absence of the superfiow {V — 0) depends only on 
|q| = 9 and has the form Sq ~ ciq + csq'^ + 0{q^) with 
ci > near q = Q. In the SMA, £q is related to the 
dynamical structure factor S{q,uj) for y = via [28t 

e„-^i-, Siq) ^ [ dujSiq,uj), (18) 



follows. When 



Vqx). 



(19) 



2S{qy 

from which S{q,uj) ^ {q^ / {2eq))5{uj — e 
F ^ 0, S'(q, uj) is given by 

S{ci,uj)^(q^/{2eq))5(uj^eq 

In spatially uniform systems, the spectral function 
/(r, uj) = I{uj) is related to the dynamical structure fac- 
tor via I{uj) = (27r)-''/dqS'(q,a;). From this and PB . 
we derive I{uj) at low uj in the following two cases, 
(i) Landau phonon instability. When the Landau insta- 
bility occurs near q = Q for C3 > 0, the Landau critical 
velocity Vc is given by the sound velocity ci. I{uj) be- 
haves as 



I{uj) oc 



UJ 



for UJ <^uj* , 



(20) 



withLj* ~ (ci-T/)3/2(c3)-i/2 andwi - {ci-V)c\'^ c^^'^ . 
(ii) Landau roton instability. We consider a system where 
the single mode Sq has a phonon-roton spectrum in the 
absence of V . The critical velocity Vc is determined by 
demanding that the condition Eq — Vcq — Q he met for 
some wave vector q = qc which lies antiparallel to the 
superfluid velocity. Around q ~ qc, Eq — Vcq can be 
expanded as ~ qc[Vc ~V) A(q — qc)^ with a positive 
coefficient A. Near V = Vc and low oj, S{q,uj) is given 
by 



Siq,uj) ^ (g2/(2£,J)^(c^-A(q-qc)- 
+ {q/{2ci))S{uj ^ ciq + Vq.^) 



qdVc 



V)) 
(21) 



4 



In the right-hand side, the first and second terms repre- 
sent roton and phonon contributions, respectively. As a 
result, I{uj) at low w is given by 



(22) 



where A and B are constants. The symbol 9 denotes the 
Heaviside step function. 

Both ([201) and show that the DBFs are enhanced 
at low bj when the superfluid velocity approaches the Lan- 
dau critical velocity. 

Stability criterion of superfluids. The enhancement of 
the DDF near the critical velocity is a common as- 
pect between a saddle-node bifurcation and the Landau 
phonon/roton instability. For those kinds of instability, 
the DDF of superfluids in d-dimensional BECs is explic- 
itly summarized in the following unified way: 



Iir,uj)^uj"fiv)+rir,uj), 
lim w~'3/'(r,w) = 

where the exponent (3 satisfies 



for ^r. 



13 = 



/3,i< d) V = V,., 
d V <Vc. 



(23) 



(24) 



The function /(r) and the exponent /3c at y = depend 
on the mechanism for destabilizing a superfluid state. 

The relation ((23|) can be also expressed in terms of 
the autocorrelation function of the local density. We 
introduce "the coarse-grained local density" n{r,t) = 
J Wa{\r — r' I )-!/)'!' (r', i)-)/)(r', i)dr', in order to avoid the 
unessential divergence of the autocorrelation at the short- 
time limit. The weight function Wa{r) (with a coarse- 
graining radius a) satisfies Wair) ^ for r ^ a(> 
0), and / Wa(|r|)dr = 1. The long-time behavior of 
the autocorrelation function C(r,t) = {h{r,t)h{r,0) + 
n(r, 0)n(r, t))/2 — {h{r,0))'^ is determined by the low lj 
property of /(r,aj) and it follows that for large t, 



C(r,i)=t-(^+i)/(r) + C'(r,t), 
\imt^^C'{r,t)t^+^ =0 



for ^r. 



(25) 



with ([Ml), /(r) is a function of r. The relation ((25)) 
reveals another aspect of instability of a superfluid state; 
a qualitative change in density correlation in the long- 
time domain is a sign that the system has reached the 
critical velocity. 

We propose ([23t and (j25l) as new criteria of stability of 
superfluids of BECs. The autocorrelation of local density 
would be more accessible than the spectral function in 
experiments of cold atoms. We thus expect that (P5|) is 
more practically useful than (P^. 

Conclusion. We examined the dynamical density fluc- 
tuations on superfluids both in the presence of the po- 
tential barrier and spatially uniform systems. On the 



basis of results for the saddle-node bifurcation and for 
the Landau instability, we propose a dynamical criterion 
for stability of superfluids. 
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